Anderson Localization of Bogolyubov Quasiparticles in Interacting Bose-Einstein Condensates 
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We study the Anderson localization of Bogolyubov quasiparticles in an interacting Bose-Einstein condensate 
(with healing length ^) subjected to a random potential (with finite correlation length ctr). We derive analytically 
the Lyapunov exponent as a function of the quasiparticle momentum k and we study the localization maximum 
fcmax- For ID speckle potentials, we find that fcmax oc 1/^ when ^ 3> ctr while fc^ax oc 1/(Tr when ^ <C ctr, and 
that the localization is strongest when ^ ~ ctr. Numerical calculations support our analysis and our estimates 
indicate that the localization of the Bogolyubov quasiparticles is accessible in experiments with ultracold atoms. 
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An important issue in mesoscopic physics concerns the ef- 
fects of disorder in systems where both quantum interference 
and particle-particle interactions play crucial roles. Multiple 
scattering of non-interacting quantum particles from a random 
potential leads to strong Anderson localization (AL) jit], char- 
acterized by an exponential decay of the quantum states over 
a typical distance, the localization length. AL occurs for ar- 
bitrarily weak disorder in ID and 2D, and for strong-enough 
disorder in 3D 13]. The problem is more involved in the pres- 
ence of interactions. Strong disorder in repulsively interacting 
Bose gases induces novel insulating quantum states, such as 
the Bose ^ and Lifshits ]|4[] glasses. For moderate disorder 
and interactions, the system forms a Bose-Einstein conden- 
sate (BEC) ]i4i,|5|,i^, where the disorder induces a depletion 
of the condensed and superfluid fractions ]0] and the shift and 
damping of sound waves 

These studies have direct applications to experiments on 
liquid ^He in porous media ]|9(], in particular as regards the 
understanding of the absence of superfluidity. Moreover, the 
realization of disordered gaseous BECs llOl 
has renewed the issue due to an unprecedented control of the 
experimental parameters. Using optical speckle fields ifTsll . 
for instance, one can control the amplitude and design the cor- 



relation function of the random potential almost at will 1 14], 
opening possibilities for experimental studies of AL iflil 17]. 
Earlier studies related to localization in the context of ultra- 
cold atoms include dynamical localization in (5-kicked rotors 
I ]_8t] and spatial diffusion of laser-cooled atoms in speckle po- 
tentials ifigli . 



Transport processes in repulsively interacting BECs can 
exhibit AL 1 1/7, However, for BECs at equilibrium, 

interaction-induced delocalizing effects dominate disorder- 
induced localization, except for very weak interactions ]|5ll3]. 
The ground state of an interacting BEC at equilibrium is thus 
extended. Beyond, one may wonder how the many -body (col- 
lective) excitations of the BEC behave in weak disorder In 
dilute BECs, these excitations correspond to quasiparticles 
(particle-hole pairs) described by the Bogolyubov theory 12111 . 
In this case, the interplay of interactions and disorder is subtle 
and strong arguments indicate that the Bogolyubov quasipar- 
ticles (BQP) experience a random potential screened by the 



BEC density ]|5|]. This problem has been addressed in the ide- 
alized case of uncorrected disorder (random potentials with a 
delta correlation function) in Ref . 12211 . 

In this Letter, we present a general quantitative treatment of 
the localization of the BQPs in an interacting BEC with heal- 
ing length ^ in a weak random potential with arbitrary cor- 
relation length (Tr. For weak disorder, we introduce a trans- 
formation that maps rigorously the many-body Bogolyubov 
equations onto the Schrodinger equation for a non-interacting 
particle in a screened random potential, which we derive an- 
alytically. We calculate the Lyapunov exponent Tk (inverse 
localization length) as a function of the BQP wavenumber k 
for ID speckle potentials. For a given ratio ^/ctr, we deter- 
mine the wavenumber fc^ax for which Tk is maximum: We 
find fcmax 1/C for ^ > ctr, and fc^ax ~ I/ctr for ^ < ctr. 
The absolute maximum appears for ^ ~ ctr, so that the finite- 
range correlations of the disorder need to be taken into ac- 
count. Numerical calculations support our analysis. Finally, 
possibilities to observe the AL of BQPs in BECs placed in 
speckle potentials are discussed. 

We consider a d-dimensional Bose gas in a potential ^(r) 
with weak repulsive short-range atom-atom interactions, char- 
acterized by the coupling constant g. Its physics is governed 
by the many-body Hamiltonian 



(1) 



+V{r)fi + {g/2)h^ - iih] 



where m is the atomic mass, ^ is the chemical potential, and 
9 and n are the phase and density operators, which obey the 
commutation relation [7i(r), ^(r')] = \5{y — r'). Accord- 



ing to the Bogolyubov-Popov theory 12 IL |23, 24], for small 
phase gradients {'t{^\V9\'^ /2m <C ^J) and small density fluc- 
tuations {5h <C ric, where = (n) and 5n ^ n — ric), 
Hamiltonian ([T]i can be diagonalized up to second order as 
H = Eq + ^t^i^' where bi, is the annihilation operator 

of the excitation (BQP) of energy e^. The many -body ground 
state of the Bose gas is a BEC with a uniform phase and a 
density governed by the Gross-Pitaevskii equation (GPE): 



H = -ft^V^(V7^)/2TO0i; + V{y) + gnc(r). 



(2) 



2 



Expanding the density and phase in the basis of the exci- 
tations, 9{r) = [-i/2VnrW]E.[/.+ (r) k - H.c] and 
5fi{r) ~ ^/nJr) X^j^l/i^C'") + H.c], the Hamihonian re- 
duces to the above diagonal form provided that obey the 
Bogolyubov-de Gennes equations (BdGE) 12511 : 



[-(nV2m)V2 + V + 2,gn, ~ /- = 



(3) 
(4) 



with the normahzation / dr [/+/^* + /^/j;*] = 25uy. 
Equations ©-(HI form a complete set to calculate the ground 
state (BEC) and excitations (BQPs) of the Bose gas, from 
which one can compute all properties of finite temperature or 
time-dependent BECs. 

Here, we analyze the properties of the BQPs in the pres- 
ence of weak disorder. According to Eqs. Q and (01), they 
are determined by the interplay of the disorder V and the 
BEC density background n^. Let ^(r) be a weak random 
potential {V ^ /i, see below) with a vanishing average and a 
finite-range correlation function, C(r) = V^c{Y/a^), where 

= ly"^) is the standard deviation, and CTr the correlation 
length of V. As shown in Refs. flSS], the BEC density pro- 
file is extended for strong-enough repulsive interactions (i.e. 
for <C -L, where £, = h/ \/4m/i is the healing length and L 
the size of the BEC). More precisely, up to first order inV^/ii, 
the GPE yields 



"-c(r) = [^J.- V(r)]/g 



(5) 



where V{r) ^ J dr'G^{r - r')Vl(r')and G^, the Green 
function of the linearized GPE IJ, |6|, r eads G^(q) = 
(27r)-'^/V[l + (|q|0^] in Fourier spacelM. Then, 



nq) = nq)/[i + (|qlO']- 



(6) 



Thus, ^ is a threshold in the response of the density Uc to 
the potential V, as V{q) ~ ^(q) for |q| < while 
V{q) < V{q) for ]q| > The potential V{r) is a 

smoothed potential IQ). If y is a homogeneous random po- 
tential, so is V, and according to Eq. (|5]l, the BEC density 
profile 71c is random but extended fl, @] . 

Solving the BdGEs (O and (01 is difficult in general be- 
cause they are strongly coupled. Yet, we show that for a 
weak (possibly random) potential V{r) this hurdle can be 
overcome by using appropriate linear combinations of the 



functions, namely g^ = ±p^^^^f+ 



with 



Pfc = a/1 + l/(fc^)2 and k = |k|. For V ^ 0, the equations 
for g^ are uncoupled [see Eqs. (|7]i and ([S])] and we recover 
the usual plane-wave solutions of wave vector k and energy 
Efe = Pfc(?i^fc^/2m) ll2lll . For weak but finite V, inserting 



Eq. (O into Eqs. (O and (0]), we find 
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(7) 



(8) 



Equations ^ and (O, which are coupled at most by a term of 
theorderof y [since |y| < \V\ and 2pk/il + pi) < 1], allow 
for perturbative approaches. Note that the functions g^ and 
g^ have very different behaviors owing to the signs in the left- 
hand-side terms in Eqs. (|7]) and Equation (O can be solved 
to the lowest order in V^/p with a Green kernel defined in 
Ref. {^]: wefindffk (r) ^ dr'G^,{r-r')V{r')g+{r'), 

where = + (fc^)^- Equation (|7]i cannot be solved 

using the same method because the perturbation series di- 
verges. Nevertheless, from the solution for g^, we find that 

1.9k /.9k I ^ li^l^l < \^\/P < 1- The coupling term in 
Eq. (|7| can thus be neglected to first order inVe^/p and we are 
left with the closed equation 



{h'/2m)V'g++Vk{r)g+ 



{h^e/2m)g+, 



where 



Vfc(r) = V{r) - 



1 



(9) 



(10) 



i + 2{kiy 

Equation (|9|l is formally equivalent to a Schrodinger equation 
for non-interacting bare particles with energy fi^k^ /2m, in a 
random potential Vfc(r). This mapping allows us to find the 
localization properties of the BQPs using standard methods 
for bare particles in ID, 2D or 3D ll27ll . However, since Vfc(r) 
depends on the wave vector k itself, the localization of the 
BQPs is dramatically different from that of bare particles as 
discussed below. 

In the remainder of the Letter, we restrict ourselves to 
the ID case, for simplicity, but also because AL is expected 
to be stronger in lower dimensions |2i|. The Lyapunov ex- 
ponent Ffe is a self-averaging quantity in infinite ID sys- 
tems, which can be computed in the Born approximation us- 
ing the phase formalism Il27ll (see also Ref. 11711 ). We get 
Ffe = (V^/8)(2m/fi2/c)2Cj.(2A:), where Cfe(q) is the Fourier 
transform of the correlation function of Vfc(z), provided that 
Ffc < fc illlBIIl. Since Cfe(g) cx {\Vk{q)\^), the compo- 
nent of Vfc relevant for the calculation of F^ is Vk{2k). From 
Eqs. dSll and ([TOll, we find 



Vk{2k)^S{kOVi2k); S{kO 



l + 2{kS,f 



and the Lyapunov exponent of the BQP reads 

Ffc = [S{kO?^k 



(11) 



(12) 
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where 7fe = {^/2^ /Z2){VJ ^lf{(T^/k'^i^)c{2ka^) is the Lya- 
punov exp onent for a bare particle with the same wavenumber 
k Hill- 

Let us summarize the vaHdity conditions of the perturba- 
tive approach presented here. It requires: (i) the smoothing 
solution (|5]l to be valid (i.e. K <C /i); (ii) the coupling term 
proportional to in Eq. ^ to be negligible (which is valid if 
K <C /i); and (iii) the phase formalism to be applicable. The 
latter requires Tk < fc, i.e. (K/m)(ctr/0^/^ < + 
1/2(A:C)2], which is valid for any A: if (K/Ai)(crR/0^^^ < 1- 

Applying Eq. (fTSl i to uncorrected potentials [C(z) = 
2D5{z) with ctr ^ 0, K ^ oo and 2D = I^^CTr / dxc(x) = 
csf\, one recovers the formula for Tu found in Ref. 12211 . 
Our approach generalizes this result to potentials with finite- 
range correlations, which proves useful since uncorrected 
random potentials are usually crude approximations of real- 
istic disorder, for which ctr can be significantly large. We 
show below that if £ < Cr, as e.g. in the experiments of 

Refs. iiS[Til[il[il Il4ll . the behavior of Ffe versus k is dra- 
matically affected by the finite-range correlations of the dis- 
order. 

Let us discuss the physical content of Eqs. (|9]l and ( fTol ). 
According to Eqs. (O and (IDi, the properties of the BQPs are 
determined by both the bare random potential V and the BEC 
density in a non-trivial way. Equation ^ makes their roles 
more transparent. As the occurrence of the smoothed potential 
V{z) in Eq. ( fTol ) is reminiscent of the presence of the mean- 
field interaction gn^^ in the BdGEs ^ and (01), it appears that 
the random potential Vfe(z) results from the screening of the 
random potential V{z) by the BEC density background 
More precisely, the expression (fTTl i for the Fourier compo- 
nent Vfc(2fc), relevant for the Lyapunov exponent of a BQP, 
shows that the screening strength depends on the wavenum- 
ber k. In free-particle regime {k ^ 1/0' we find that 
the Lyapunov exponent of a BQP equals that of a bare particle 
with the same wavenumber {Tk — Jk), as expected. In the 
phonon regime (k ^ 1/0' the disorder is strongly screened 
and we find ^ 7^, as in models of elastic media 12811 . 
Here, the localization of a BQP is strongly suppressed by the 
repulsive atom-atom interactions, as compared to a bare par- 
ticle in the same bare potential. These findings agree with 
and generalize the results obtained from the transfer matrix 
method, which applies to potentials made of a ID random se- 
ries of (5-scatterers ll22ll . 

Our approach applies to any weak random potential with a 
finite correlation length. We now further examine the case of 
ID speckle potentials used in quantum gases ifloifTl . 12. 13. 
[Till . Inserting the corresponding reduced correlation function, 
c(k) = \/W2(l - k/2)6(1 - k/2) where 9 is the Heaviside 
function llTtj, into Eq. ( fT2] ), we find 



a,k^l-ka,) 



■ my 



(13) 



M/ [1 

which is plotted in Fig. [T] To test our general approach 
on the basis of this example, we have performed numeri- 
cal calculations using a direct integration of the BdGEs ^ 
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Figure 1: (color online) Density plot of the Lyapunov exponent of 
the BQPs for a ID speckle potential. 



and (HI in a finite but large box of size L. The Lyapunov 
exponents are extracted from the asymptotic behavior of 
\og[ri.{z) /rk{zk)]/\z — Zk\, where z^ is the localization center 
and 7'A; (z) is the envelope of the function g'^, obtained numer- 
ically. The numerical data, averaged over 40 realizations of 
the disorder, are in excellent agreement with formula (T3[ as 
shown in Fig.|2] These results validate our approach. It should 
be noted, however, that our numerical calculations return BQP 
wave functions that can be strongly localized for very small 
momenta k. This will be discussed in more details in a future 
publication 129]. 

Of special interest are the maxima of Fj,, which denote a 
maximum localization of the BQPs. It is straightforward to 
show that, for a fixed set of parameters (Vj,/^, ctr), F^. is 
non-monotonic and has a single maximum, fctnax, in the range 
[0, I/ctr] (see Fig.|2l). This contrasts with the case of bare par- 
ticles, for which the Lyapunov exponent -jk decreases mono- 
tonically as a function of fc, provided that c(2/s(Tr) decreases 
versus k (which is valid for a broad class of random poten- 
tials 12711 ). The existence of a localization maximum with re- 
spect to the wavenumber k is thus specific to the BQPs and re- 
sults from the strong screening of the disorder in the phonon 
regime. In general, the value of k^.^^, plotted in the inset of 
Fig. Inversus the coiTelation length of the disorder, depends on 

both ^ and ctr. For CTr < we find /cmax ^ ^ (l - , 
so that the localization is maximum near the crossover be- 
tween the phonon and the free-particle regimes as for un- 
corrected potentials ll22ll . For (Jf^ ^ ^ however, we find 
fcmax — 2/3(Tr, so that fcmax is no longer determined by the 
healing length but rather by the correlation length of the dis- 
order, and lies deep in the phonon regime. For k > I/ctr, F^ 
vanishes. This defines an effective mobility edge due to long- 
range correlations in speckle potentials, as for bare particles 



Finally, let us determine the absolute localization maxi- 
mum. The Lyapunov exponent Ffe decreases monotonically 
versus ^ and K/m- However, for fixed values of Vr/p and 
F/j has a maximum at ctr = ■\/3/2 ^ and k — ^~^/-\/6 
(see Fig.[T]i and we find the corresponding localization length 
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Figure 2: (color online) Lyapunov exponent of the BQPs for a ID 
speckle potential. The lines correspond to Eq. l |13l l and the points to 
nu meri cal results for ^ = 1.25 x 10"'' L, Vr = 0.075/i and ctr — 
•\/3/25 (solid red line), (Jr — 3.7^ (dashed blue line), ctr = 0.4^ 
(dash-dotted green line). Inset: localization maximum versus the 
ratio of the correlation length of the disorder to the healing length of 
the BEC (the dash-dotted and dashed lines correspond to the limits 
Ur <^ ^ and ffR ^ 5 respectively; see text). 



(in 



l/rniax)- 



,(0 = (512V6/9^)(/./K)'e 



(14) 



At the localization maximum, we have CTr ~ ^, so that the 
disorder cannot be modeled by an uncorrelated potential, and 
the long-range correlations must be accounted for as in our 
approach. For ctr = 0.3/xm 1114(1 and Vr = 0.2fj,, we find 
in,ax — 280fim, which can be smaller than the system size in 
disordered, ultracold gases 

In conclusion, we have presented a general treatment for the 
AL of BQPs in an interacting BEC subjected to a random po- 
tential with finite-range correlations. We have calculated the 
Lyapunov exponents for a ID speckle potential and we have 
shown that the locaUzation is strongest when CTr ~ ^. We 
have found that the localization length can be smaller than the 
size of the BEC for experimentally accessible parameters. We 
expect that the AL of BQPs could be observed directly, for 
instance as a broadening of the resonance lines in Brag g sp ec- 
troscopy, a well mastered technique in gaseous BECs [320. 
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